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Abstract. This article reports on a student summer project performed in 2006 at 
the University of Frankfurt. It is addressed to undergraduate students familiar with 
the basic principles of relativistic quantum mechanics and general relativity. The 
aim of the project was to study the Dirac equation in curved space time. To obtain 
the general relativistic Dirac equation we use the formulation of gravity as a gauge 
theory in the first part. After these general considerations we restrict the further 
discussion to the special case of the Schwarzschild metric. This setting corresponds to 
the hydrogen atom, with the electromagnetic field replaced by gravity. Although there 
is a singularity at the event horizon it turns out that a regular solution of the time 
independent Dirac equation exists. Finally the Dirac equation is solved numerically 
using suitable boundary conditions. 
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1. Introduction 



The most important success in the early days of quantum mechanics was the appropriate 
description of the hydrogen atom. It is not exaggerated to say that quantum mechanics 
was discovered in its final form in the 1920ies in the attempt to describe the hydrogen 
atom which could not be understood in terms of a classical theory. 

To describe the behaviour of an electron with the help of quantum mechanics means 
to solve the time independent Schrodinger equation (respectively the Dirac equation in 
the relativistic case) with a coupled Coulomb potential and thus to find the eigenstates 
of the corresponding Hamiltonian. The existence of bound states within an atom 
has fundamental importance since it is the reason for the stability of matter. So 
it is suggestive to consider a similar situation, namely the behaviour of an electron 
within a spherical symmetric gravitational field. Since gravitation is much weaker than 
electromagnetism (about a factor of 10 43 in case of the mass and the charge of an 
electron) the body being the source of the gravitational field must have a huge mass to 
lead to an effect of the same order of magnitude as in the electromagnetic case. The 
most massive objects known today are black holes. 

Although the situation is in a certain sense similar to the case of the hydrogen 
atom, to describe gravitation one has to use general relativity instead of Maxwell 
theory. Therefore, the problem becomes much more involved because of the fundamental 
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difference between the framework of electromagnetism and general relativity. In the case 
of general relativity one cannot simply couple a potential to the Hamiltonian but the 
structure of space time itself is changed by the presence of a massive body. This leads to 
a metric tensor g^ u (x) which differs from the usual metric of flat Minkowski space time 
^.v _ diag(+l, — 1, — 1, — 1). The structure of space time near a spherical symmetric 
body is described by the Schwarzschild metric. In the case of the Klein-Gordon equation 
the metric appears explicitly in the equation and one has to replace the Minkowski 
metric by the Schwarzschild metric to obtain the correct equation adapted to the new 
situation. The situation is very different in the case of the Dirac equation, because 
the metric tensor does not appear explicitly in the Dirac equation. Besides this the 
Dirac equation does not describe a scalar valued state but a Spinor. I.e. that we have 
to incorporate the effect of the new metric on the Spinor space. So the situation is 
more difficult than in the case of the Klein-Gordon equation where the new metric has 
just to be inserted. Therefore we will have to derive the general form of the Dirac 
equation in curved space time by considering the gauge formulation of general relativity 
as mentioned above. This is the subject of the second section. After this we will explore 
the special case of the Schwarzschild metric and look for a suitable representation to get 
an equation which is as simple as possible. 

However the changing of the metric must have an effect similar to the coupling 
of a potential. Indeed, as will be shown below, the Dirac equation takes the form of 
the free Dirac equation with an additional perturbative term in the Hamiltonian when 
a certain representation is chosen. It remains the task to solve the equation. The 
main reason why it is more difficult to solve this problem compared to the hydrogen 
atom is the singularity at the event horizon in addition to the singularity at the origin. 
However, using a series expansion it turns out that a solution can be found which is 
regular at the horizon. Further it can be shown that a certain boundary condition has 
to be fulfilled there. Although the probability density of the particle does not vanish 
within the event horizon we will omit this area. With the help of the above mentioned 
series expansion and boundary condition around the horizon we will get an initial value 
infinitesimal close to the event horizon from where a numerical integration procedure 
can be started. To start this integration procedure not directly at the event horizon 
is necessary because of the possibility to get to a solution which is not regular there. 
As in case of ordinary quantum mechanics the solution has to approach zero at infinity 
to represent a regular bound state that is normalizable. With this condition we can 
find the eigenstates by varying the energy as long as the wave function does vanish at 
infinity. The stationary states have complex energy eigenvalues as we are dealing with 
a non-hermitian Hamiltonian. The complexity of the energy eigenvalues is connected 
to the instability of the states and with the help of the imaginary part of the energy 
it is possible to determine the decay times. However, an exact prediction of the decay 
times exceeds the scope of this article. A discussion of this topic can be found in ref. 
p. The topic of the summer project as presented here is based on refs. [H El Ej and 
therefore owes many ideas to these papers where the present problems have been treated 
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for the first time. A similar situation concerning the Klein-Gordon equation has been 
considered in ref. [I]. 

2. Implementation of general relativity into the Dirac equation 

Since we are going to study a fermion in a gravitational field we have to consider 
the coupling of the Dirac equation to gravity. According to the other fundamental 
interactions occurring in nature gravity can also be formulated as a gauge theory. (For 
an introduction to the idea of gauge theories in general we suggest ref. [5] (Chapter 15).) 
The difference is that not local invariance under internal symmetric transformations is 
demanded but local invariance under Lorentz transformations is considered. 

2.1. Notation 

The following symbols will be used: 

ip = Dirac Spinor (in coordinate representation). 

= vector in Hilbert space describing an arbitrary quantum state and = the 
corresponding dual vector (refraining from any special representation). 

d = -2- 

ryt* = Dirac matrices. 

= Lorentz transformation matrices and (A -1 )^ = the corresponding inverse 
transformation. 

U = arbitrary unitary operator and U' = its hermitean conjugated operator. 

Besides we will make the arrangement that Greek indices refer to arbitrary global 
coordinates and Latin indices refer to the local coordinate system in which the metric 
tensor takes the form of the special metric tensor of flat Minkowski space time r) mn . 
If a certain index adopts a specific value the special local coordinates are denoted with 
a bar to distinguish them from arbitrary global coordinates. As usual natural units will 
be used which means that h — c — G — 1. 

2.2. The general shape of the equation 

Like in case of usual gauge theories one starts with the Lagrangian of the free field 
equation, in our case the Dirac equation (first published in ref. [6]) 

C = tfli'fdr - m)ip. (1) 

In the framework of relativistic quantum mechanics a Lorentz transformation is 
represented by an unitary operator which acts on a state in Hilbert space. Infinitesimal 
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transformations are provided by the generators of the group. The generators of the 
Lorentz group obey the following Lie Algebra 

Pa/3, T, y s\- = Tj^pLaS — T] 7a Eps + TjspTjSa ~ (2) 

If we want to analyse the action of Lorentz transformations on quantum mechanical 
states of fermions described by the Dirac equation we have to consider the Dirac Spinors 
representation. In this representation the generators are given by 

£^ = -[7^7*]- ( 3 ) 

(For an extensive treatment of the theory of group representations in the context of 
physics including a derivation of Eq. (j2J) and the Dirac Spinor representation of the 
Lorentz group we suggest ref. [7] (sections 2.2-2.4 and section 5.4, resp.).) 

It is a priori clear that the Dirac equation is invariant under arbitrary global Lorentz 
transformations, i.e. independent of the space time point because it was derived from 
this constraint. (An explicit proof of this invariance property can be found in ref. [8].) 
In contrast to this, the well known Dirac Lagrangian is not invariant under local Lorentz 
transformations, but these are exactly the ones needed in the case of general relativity. 
In general relativity flat Minkowski space time is replaced by a curved space time (this 
topic is originally presented in ref. [9]). This means that at every point of space time the 
metric tensor g^ v looks different. If one wants to compare two mathematical objects at 
different space time points, one first has to transform to the frame in which the metric 
tensor looks like the metric tensor i] fiu of flat Minkowski space time. Thus a specific 
Lorentz transformation at each space time point has to be performed. This is the reason 
why we have to consider local Lorentz transformations within the framework of general 
relativity. Thus, it is our next task to formulate the Dirac equation in a way that obeys 
local symmetry transformations 

i) -> U(x)ij , «9 M -> A;(x)d„. (4) 

A Lagrangian of the following shape satisfies invariance under local Lorentz 
transformations 

£ = i){%YD^ - m)if>, (5) 

where denotes a new derivative operator (covariant derivative) which transforms 
under local Lorentz transformations according to 

Dfj, -> N-;U{x)D v U\x). (6) 

It is clear that such a covariant derivative obeys local Lorentz invariance, because the 
U operators to the left and right of the derivative operator cancel the operators which 
act on the right and on the left state in the Lagrangian. 

The covariant derivative applied to the metric tensor is zero. This means that 
is the derivative connecting different space time points endowed with a local pseudo- 
Euclidean coordinate system in which g^ u = r/^. Since we are regarding Spinors on the 
space time it is also necessary to consider the action of on the Dirac Spinor space. 
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Now we have to find out how this covariant derivative looks like. To do this, we 
assume that is of the following shape 

D m = e^d„+iu li ). (7) 

Here we have introduced the vierbein e v m describing the transformation between arbitrary 
global coordinates and the local coordinates in which the metric tensor looks like the 
metric tensor of flat Minkowski space time. This means that the following relation is 
valid 

Vmn e m e nfi , M^' (^) 

which can be seen as the definition of the vierbein. It expresses that if we know 
the vierbein at each space time point we can obtain the metric field meaning that 
it contains the same information about the structure of space time. So the vierbein 
represents another way to describe the gravitational field and can replace the metric 
tensor. Further, the vierbein fulfils the relations 

(An extensive introduction to the vierbein-formalism can be found in ref. |10j). 

By inserting Eq. (jZ|) into the above transformation condition Eq. ([6]) and rewriting 
it a little bit one gets 

+ zav) -> A£e£ [8, + U(x)d,U\x) + iU(x)u^(x)} . (10) 

So the covariant derivative e^c^ + iu^) obeys condition (J7J) if and uj^ transform 
according to 

e»^K& , u^Uu^-iUid^). (11) 

The matrices describe the Dirac Spinor connection. They can be performed by the 
above generators of the Lorentz group of Eq. (E]) 

w„(x) = ±u™(x)Z mn . (12) 
So we have found the Lagrangian which is invariant under local Lorentz transformations 



£ = 4> 

leading to the Dirac equation in general relativity 



V, (13) 



ip = 0. (14) 



(For an extensive discussion of gravity as a gauge theory the reader is referred to ref. 

mao 

Combining the well-known anti-commutation relations for the Dirac-matrices 
[7^,7^] + = 2t]^ u and (jHJ) leads to the new anti-commutation relations 

K 7 m ,<7 B ]+ = 2<r, (15) 

for the Dirac matrices in curved space time. We will need this relation in the next 
chapter. 
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2. 3. Calculation of the Spinor coefficients 

Now we have to determine the specific form of the Spinor connection coefficients uj™ n . 
It is important to notice that our coefficients cj™ n contain one index (//) referring to 
the global coordinate system and two indices (m, n) referring to the local coordinate 
system where g mn = r\ mn . First, we express the above coefficients u™ n by the usual 
Christoffel symbols Y p defining the parallel transport between vectors in the tangential 
vector space of different space time points. They are defined by the condition that the 
metric shall be left invariant by the covariant derivative 

D»(g P a) = (16) 

and have the following shape 

= \g p °{d v g^ + d,g ua - d a9txv ). (17) 

(The explicit derivation of this relation can be found in any book on general 
relativity, e.g. [TT].) 

We can express the covariant derivative of the vierbein by applying the Christoffel 
symbol to each index leading to 

D e m = d e m + T m e n - Y p e m (18) 

Inserting Eq. ([8]) into Eq. (|T6|) and considering the case v = p one obtains 



D,{g uu ) = D^e n v r, mn ) = D^)\ mm ) 

= VmmD^e™) 2 = 27? mm e™ZVr = 0. (19) 

which means that the covariant derivative applied to the vierbein gives zero 

= 0. (20) 

Here we have used that D^rjmn = and i] mn is a diagonal matrix. Combining Eq. (1181) 
and (120]) results in 

d,e™ + r™x - r>™ = o. (21) 

By moving the second term of Eq. (12 ip to the r.h.s. side of the equation, multiplying 
with e v n and raising the flat coordinate indices with 7] al3 we arrive at a relation between 
the above Dirac Spinor coefficients and the Christoffel symbols in arbitrary coordinates. 

< n = IT = -V np e;(d,e™ - I> p m ). (22) 

Here, we can identify the coefficients of the connection uj™ n with the Christoffel 
symbols expressed in terms of the local coordinates, because the connection between the 
tangential vector spaces of different points on the space time manifold is independent 
of the special space by which it is represented. If one recalls now that the vierbein 
describes the transformation between global and local coordinates, one can express the 
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coefficients by the vierbein. Using Eq. (jS]) in Eq. (fT71) . inserting it to Eq. (T2"2"j) and 
transforming the obtained expression a little bit, we arrive at 

u rnn = 2e ^d^ + e^e^ . (23) 

So we have expressed the connection coefficients in terms of the vierbein. (The derivation 
of this formula can also be found in [12]. For advanced aspects concerning the relation 
between the usual formulation of general relativity and the vierbein formalism we suggest 
Refs. [ID], [12] and [13].) 



3. Solution of the Dirac equation in a Schwarzschild background 



3.1. The choice of a suitable formulation of the Dirac equation in a Schwarzschild 
background 

So far we have formulated the Dirac equation in a general space time with an arbitrary 
metric tensor and an arbitrary connection between different space time points. To be 
specific we will from now on focus on the special case of the Schwarzschild metric, 
because we are investigating a particle in a spherically symmetric gravitational field. 



The Schwarzschild metric (first presented in 
following way 

( 

' TT^MT 






is usually represented in the 



V 



















(24) 



— r 2 sm(8) J 



A spherically symmetric object only influences the radial part of the metric. Therefore, 
we have to deal with the time and the radial part of the metric and do not consider the 
angular coordinate. This leads to 



e 5 7 5 + e° 7 f 
e r 5 7° + e r ,Y 



(25) 



The components of the Schwarzschild metric are not fully determined due to the fact 
that the metric can be expressed as a function of a new time coordinate according to 
a transformation of the form t — > t + x(r) where x denotes an arbitrary function of r. 
It is important to note that this transformation of the time coordinate indeed changes 
the appearance of the metric, but does not affect the goo-component of the metric. This 
component is equal in all representations 



#00 



(26) 



Besides this, the vierbein is still under-determined if a special time coordinate is chosen 
because of the possibility to apply a radial dependant Lorentz boost to the system 
connected to the choice of a global gauge. 
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This means that two degrees of freedom remain. However, two constraints can be 
formulated. The first constraint can be derived from the fact that the new objects e° 7 M 
and e^7 M have to fulfil the same Lie- Algebra as 7 and Y which means that 

[ 7 °,7 r ]- = [e^7 m ,e;7 n ]- (27) 

And the second constraint arises from the invariance of the goo-component of the metric. 
By applying the covariant metric tensor g^ v to Eq. (IT51) we get two further equations 

K7 m ,<7„] + = 2#, 

K7m,e n uln } + = 2 9ilv . 1 ] 

So we can formulate the two following constraints 

pOpT _ pOpT 1 

[el? -{elf = (1-2M). ^ 

It is very important to emphasise that all choices of the vierbein fulfilling 
the constraints ( 1291) represent the usual Schwarzschild metric but in different 
representations. It turns out that the Dirac equation in a Schwarzschild background 
becomes quite simple if we choose the coefficients as 



eg = l » 4 = , 4 = 1 , el = -(^y. (30) 

It is left to the reader to show that these choices fulfil the conditions (1291) . The inverse 
components of the vierbein are determined by Eq. (jHJ) leading to 



/2M\5 
V r 



< = 1 , e° = 0. (31) 



By the special choice of the vierbein in Eq. (|30|) . the special form of the metric is 
completely determined by Eq. ({151) and Eq. (|28|) . To find the specific Dirac equation 
we have to specify the connection coefficients io™ n by inserting the above choice of the 
vierbein to Eq. (T23]) . In the expression arising from this most terms vanish and one 
obtains the following coefficients 

The coefficients with equal upper indices vanish because of the antisymmetry. By 
inserting Eqs. (130]) and (1321) to the general Dirac Eq. (1141) and using the antisymmetry 
of the generators and the cooresponding coefficients leading to u;™ n S mn = 2o;° r £of we 
get the following expression 

(^-m-z 7 °(— J 9r-T r (-2 (—J - Ufa, 7r]-)V> = 0.(33) 
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7° and y r obey the same anticommutation relation as 7 and 7* where i runs from 1 to 
3. From this we see that [70, 7r]- = — 27- r 7o an d further that 7 r 7 r = 1. So we find that 
in this representation the Dirac equation takes the following form 

1 

(iy^dp -m-ty° * (d r + |-))V> = 0. (34) 

(This result for the Dirac equation in a Schwarzschild background with our special 
choice of gauging corresponds exactly to the one obtained in pQ). 

3.2. Reformulation to a system of differential equations 

In the previous section we have derived the Dirac equation on a Schwarzschild 
background in a suitable representation. It has the form of the free Dirac equation with 
an additive radial dependant perturbation term. Now the task remains to reformulate 
the free part of the equation and to express it in spherical coordinates and then to 
separate the radial part of the state from the spherical symmetric one. 

First we rewrite Eq. ( IMl) by expressing the 7-matrices by the well-known a- and 
/3-matrices which are defined by the relations 70 = (3 and y k = (3a k for k = 1, 2, 3 and 
obtain 

(i(3d t + t(3a k d k -m-i/3 (—) * (d r + |-))^ = 0. (35) 

\ r / 4r 

Rearranging this equation and multiplying it with (3 we get 

{ia k d k )ij = {-id t + f3m + i (—) 2 (d r + i)V> = 0. (36) 

V r / 4r 

Since we are finally interested in stable orbits (the energy eigenstates) we separate the 
time to get the time independent Dirac equation. 
Choosing the following ansatz 

^ = e - iEt V, (37) 

leads to 

1 

{ia k d k )^ = (-E + /3m + i (—) ' {d r + (38) 

V r J 4r 

We can rewrite this to 

H* = {H Dirac + H INT )V = EV, (39) 



where 



and 



H Di rac = (~ia k d k + [3m), (40) 



Hint = i 



r 
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\1/ shall describe an eigenstate to J 2 , J z and P where J describes the norm of the overall 
angular momentum which means the sum of the orbital angular momentum L and the 
Spin S, J z its component in ^-direction and P the parity. 

The Parity will be denoted by the quantum number uo (not to be confused with the 
Spinor coefficients uo™ n ). We set uo = +1 for states with parity (— ) J+ 5 and uo = — 1 for 
states with parity (— ) J ~^. 

In the following yffj will describe an eigenstate to J 2 and J z where the parity is 
given by (— ) L which means that L can adopt the two values L = J + \uo and L = J— \uo 
and that yffj and y^j have opposite parity. 

Next we make the following ansatz for \l/ 



(42) 



where $ and x describe arbitrary radial states. By following a derivation also used in 
case of the relativistic hydrogen atom (where the free part of the Dirac equation has 
to be reformulated in spherical coordinates also and which can be found in Ref. [15] 
for example) and splitting the two components of the equation into a system of two 
differential equations we obtain 



— d r 



d r + 



w{J+\) 



w{J + 



X 



.{2My l, 



, 2M Y 2 (x ^ 



X- 



(43) 



These are the differential equations for a Dirac particle in a Schwarzschild metric 
we have to solve. 



3.3. Solution of the time independent equation with suitable boundary conditions 

After deriving the system of differential equations the final task is to calculate the energy 
eigenstates by solving the set of coupled partial differential equations. Because the set of 
equations contains complex numbers it is convenient to split the two radial components 
$ and x m a real and imaginary component resulting in four equations which can be 
integrated numerically. 

However, before we do this, we have to make sure that there exists a regular solution 
at the event horizon and to determine the necessary boundary conditions. Therefore 
we will transpose the above equation system. In both Eqs. (j4"3"j) appear the derivatives 
of $ and x- By solving one equation for <9 r $ and the other for d r x and inserting each 
equation into the other we get the following system 



2M\ 



(2M\ 



+ i I— [E + m + 



2M\ 
r J 



1] 

4r J 
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2M\ 



drX 



,2_My 

V r / r 4r J 



2M\ 2 



[E 1 — m + i 



+ 



2MY . w{J+\) 1 

— ^ + —+(£ + m 



v 



2M^ 1 



(44) 



After this transformation the singularity at the event horizon (r s = 2M) shows up 
explicitly. Nevertheless it is possible to find a regular solution. To obtain this solution 
we perform a series expansion around the event horizon 



$ = X>(r-2M)* , X = X>( r - 2M )' 



(45) 



i=l i=l 

By inserting now the series expansion into Eq. (j44p we see that it is a solution of our 
differential equation system at the event horizon r = 2M under condition that 



%E - 




Here we have used 

oo 

X>(r-2M)* 



=2M 



(46) 



(47) 



and 



2M\ 



7 1=1 



r 

J oo 



r - 2M) 
- 2M) (l - 1) 
r - 2M ) ' 



2M 



2M ~ / 



=2M 



\r=2M 



0. 



(4f 



To begin an integration of the differential equation system Eq. (1441) one needs the start 
values of $ and x- By using the series expansion Eq. f|45l) we have found one constraint 
for the start values of the integration of <fi and \. The remaining degree of freedom is 
only connected to the the overall normalisation of the wave function. It is determined 
by the normalisation condition that the integral over the squared norm has to be equal 
to unity because of the probability interpretation of the wave function. So the wave is 
defined apart from a physically irrelevant overall phase. 

The integration cannot be started directly from the horizon because of the problem 
with the non-regular solutions. However, if we consider our series expansion Eq. (1431) 
to the first order infinitesimal adjacent to the event horizon we can find the coefficients 
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a\ and b\ by inserting it to the equation system f!44p and solving it. Because of the 
complexity of ai and b\ we have effectively four equations. The calculation is quite long 
but mathematically trivial. We will therefore refrain from presenting it here. 

The coefficients ai and b\ of the first order series expansion Eq. fj4"5l) together 
with the start values of <ft and \ define the wave function infinitesimally close to the 
horizon from where the numerical integration can be started. The energy eigenvalues 
are found by the condition that the wave function has to approach zero at infinity to 
assure integrability. 



3.4- Non-hermiticity of the Hamiltonian 

Unfortunately, we cannot be sure that the eigenvalues of the Hamiltonian in Eq. ( 1391 ) 
H = Hpirac + Hj nt are real. In general real eigenvalues are only obtained from 
hermitean operators. In the case discussed here, we will show that the Hamiltonian 
is not hermitian. The Hamiltonian is the generator of the group describing the way 
quantum mechanical states in Hilbert space evolve with time. This means that 

U{t 2 ,t 1 )=exp[-iH(t 2 -t 1 )}, (49) 

where U(t 2 ,ti) describes the relation between a quantum mechanical state \if)(ti)) at 
time t 1 and the same state \ip(t 2 )) at time t 2 within the Schrodinger picture, i.e. 

m 2 )) = u{t 2 M)m l )). (50) 

If the time evolution operator U(t\,t 2 ) is unitary this results in a constant overall 
probability. But Eq. P9|) tells us that this is only the case because of the hermiticity 
of the Hamiltonian. In the case of a non-hermitean Hamiltonian unitary breaks down 
and the norm of the state can change. In our case this is due to the fact that the 
singularity at the centre of the black hole leads to an instability of the states resulting 
in a probability current through the event horizon. Before drawing our attention to 
this phenomenon again we will first show this non-hermiticity of the Hamiltonian of Eq. 

Hermiticity implies that the expectation value of an operator is real 

(V\H\V) = (V\H\V)*. (51) 

In coordinate representation the inner product of Dirac Spinors is defined as 

00 00 

—00 —00 

As the free Dirac Hamiltonian Hoirac is hermitian, we have to check the additional 
perturbation Hamiltonian Hj nt 

00 1 
(*\Hmtm = J ^7°H7° (^r) 2 {Or + ^ d 3 r (53) 
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mHmtW = J V> T 7°(*7° (^) 2 (d r + d 3 r. (54) 

— oo 

If we integrate the derivative term by parts and use the condition that the norm of the 
wave function vanishes at the boundary we find the expression 

<*i*«i»>- = 4* / ((w) (-i (^)*) r + v (i i^ff^) «•) r 2 * 

i 00 i 

+ 4vr ^ T (z (2Mr 3 ) 5 j ^* |r=o +4vr ^ ^ T (2Mr 3 )'j ^* dr. 



(55) 



Here we have used that (70) 2 = 1- If we now use that = ip T <j) we obtain 



00 

+ 4vr y V>V (-i 7 ° (d r (2Mr 3 ) 5 - (2Mr 3 ) 5 — ^ V 


+ 4tt ^ fz (2Mr 3 ) 5 j ^ |r=o 

00 

= ($\H Int \^) + i J 5(r)^ (2Mr)* ip d 3 r. (56) 

—00 

So we have 7^ H and the energy eigenvalues need not to be real and so we have to 
extend our search to complex numbers. 

However if we have found the eigenstate we are mainly interested in the real part 
in the context of this article. The real part of the energy has to be smaller than the rest 
energy of the particle to provide a (if we abstain from instability) stationary state. It 
turns out that the imaginary part is always negative. Therefore we define 

E = E R - iE Im , (57) 

with Ej m > 0. As already mentioned in the introduction the imaginary part is connected 
to the instability of the states. This instability corresponds the analogue fact in the 
framework of particle mechanics that a particle which has passed the horizon cannot 
get back to the area outside the horizon and disappears in the singularity at the centre 
of the black hole. In the case of quantum mechanics this means there is an effective 
inward current. By recalling now Eqs. (115)1 and (1501) and using the eigenvalue equation 
we obtain 



\i/)(t 2 )) = exp[-iE(t 2 -ti)]\iJ){t 1 )). 



(58) 
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With the help of Eq. (1571) we can transform Eq. (I58p to 

|V(t 2 )> = {exp [-iE R (t 2 - h)] |V(ti))}exp [-E Im (t 2 - h)] . (59) 

So it can be seen directly that a negative imaginary energy part results in an exponential 
decrease of the probability density This is only another way to express the fact that 
probability density vanishes with time. Such a phenomena can only be interpreted as 
the effect of a current sink at the singularity of the centre. If probability density extends 
to the area within the horizon it cannot escape and has to end at the singularity in the 
centre. Without this special structure of space time connected to black holes such a 
strange phenomena could not take place. This means that the decay rate of the states 
depends on the properties of the wave function near the origin. A small probability 
density near the singularity leads to a long decay time. Thus states with a large 
angular momentum, i.e. a small probability density near the origin are very stable. 
This shows up explicitly in the decreasing imaginary part of the energy as a function of 
the displacement of the probability density outwards (see Figures) . For a more extensive 
discussion of this decay phenomenon the reader can consult ref . [I] . 

4. Figures and Discussion 

The solutions of the Dirac equation in a Schwarzschild background differ significant from 
the solutions of the corresponding equation with a coupled Coulomb-potential describing 
a hydrogen atom. The main difference is the non-stability of the states connected to 
the non-Hermiticity of the Hamiltonian. In the gravitational field real bound states 
exist only in an approximative sense and if the angular momentum becomes large the 
main part of the probability density lies far away from the origin. Using the Newtonian 
gravitational potential would correspond exactly to the hydrogen atom with other pre- 
factors. Thus the deviations in case of gravity arise from the description obtained by 
considering general relativity. 

The other comparison concerns the relation to the orbits of classical particles. In case 
of classical mechanics on a Schwarzschild background there are two different kinds of 
states or orbits corresponding to the obtained quantum mechanical states. On the one 
hand there are orbits with high angular momentum preserving particles from falling into 
the black hole. On the other hand there are possible orbits making the particle crossing 
the horizon. If such a crossing of the horizon takes place the particle can never reach the 
area outside again and falls to the centre of the black hole. The orbits or states of the 
first category are stable in an absolute sense. This is different in the case of quantum 
mechanics where no such discrete classification can be made leading to the fact that all 
states have an imaginary part which never vanishes. 

Qualitatively, we might also obtain similar results with the Klein-Gordon equation 
since the special properties of the considered situation, e.g. the instability of the 
states, are independent of specific aspects of different particles. However, one does 
not expect to get the same energy levels because of the special aspects related to the 



CONTENTS 



16 



inner degree of freedom of Fermions leading to an additional Spinor contribution to 
the Hamiltonian. The additional Spinor connection term gives the reason why the 
energy eigenvalues for Dirac particles are different from those of Klein-Gordon particles. 
Besides, the description of fermions on curved space time is only possible by introducing 
the mathematical concept of the vierbein. In spite of the more difficult solution it is 
more interesting to consider the Dirac equation because most fundamental particles are 
fermions and thus obey the Dirac equation. 

To obtain concrete Figures and energy eigenvalues we have chosen the special case of 
an electron having a rest mass of m e = 9.109 • 10~ 31 kg in the gravitational field of a 
black hole with mass M = 1.825 • 10 14 kg. To specify a particle with a certain rest 
mass is important since the concrete eigenvalue of the energy depends on it according 
to the relativistic relation between energy and momentum. This setting is similar to an 
electron in the electromagnetic field of a proton, hence a hydrogen atom. The mass of 
the black hole chosen here corresponds to a Schwarzschildradius of r s = 2.7 • 10~ 13 m. 
For the lowest state of J = 1/2 (see Figured]) i.e. without orbital angular momentum 
and the lowest states for J = 3/2 and J = 5/2 (see Figures [2] and [3]) |$| 2 + \ x\ 2 is plotted 
as a function of r/r s . One can see that increasing orbital angular momentum leads to 
a distinct movement of the probability density outwards. In Figure [1] the probability 
density at the horizon is still more than 50 % of the peak probability density. In the 
case with orbital angular momentum of 1 one observes that the contribution near the 
horizon nearly vanishes. 
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Figure 2. |$| 2 + \x\ 2 as function of r/r s for the lowest state of J = 3/2 (lp-state) 
corresponding to an energy of Er — 9.83 • 10 _1 m c c 2 , Ei m — 3.4 • 10~ 5 m o c 2 . 
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